Combination quantum oscillations in canonical single-band Fermi liquids 
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Chemical potential oscillations mix individual-band frequencies of the de Haas- van Alphen (dHvA) 
and Shubnikov-de Haas (SdH) magneto-oscillations in canonical low-dimensional multi-band Fermi 
liquids. We predict a similar mixing in canonical single-band Fermi liquids, which Fermi-surfaces 
have two or more extremal cross-sections. Combination harmonics are analysed using a single-band 
almost two-dimensional energy spectrum. We outline some experimental conditions allowing for 
resolution of combination harmonics. 
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Magnetic quantum oscillations of magnetisation 
(dHvA effect) and resistivity (SdH effect) are unequivo- 
cal hallmarks of the Fermi-liquid, providing most reliable 
and detailed Fermi-surfaces in particular in layered 
organic metals @, 0] and almost two-dimensional (2D) 
superconductors like Si'2Ru04 An interesting feature 
of dHvA/SdH oscillations is a measurable difference be- 
tween canonical and grand canonical ensembles, which is 
most pronounced in multi-band low-dimensional metals 
Q. Thermodynamically two ensembles must be identi- 
cal, but quantum fluctuations are fundamentally different 
depending on whether measurements are performed on 
either closed or open system with fixed electron density, 
n e , or chemical potential, /i, respectively. The difference 
between corresponding free energies is tiny, since it is 
proportional to fluctuations of the carrier density, how- 
ever, the effect on quantum corrections in magnetisation 
and conductivity is significant. 

In particular, there are combination frequencies in 
dHvA/SdH oscillations of a two-dimensional multiband 
metal with fixed n e , predicted by Alexandrov and 
Bratkovsky (AB) @, and studied numerically @, @, 0, H, 
HG3 and analytically (HI El El . The effect was exper- 
imentally observed in different low-dimensional systems 
0, El. Obviously, there are no chemical potential 
oscillations when fi is fixed by a reservoir, so there is 
no mixing of the individual-band fundamental frequen- 
cies in the Fourier transform (FT) of magnetisation in 
an open (grand-canonical) system. Importantly, samples 
are normally placed on non-conducting substrates with 
no electrodes attached, so the system is closed in actual 
dHvA experiments. 

As it happens the fundamental frequency mixing due 
to the chemical potential oscillations (AB effect) may be 
obscured by mixing due to the magnetic breakdown [l6j 
(MB effect), as discussed by Kartsovnik Q. The MB ef- 
fect is the switching of two close electron orbits in differ- 
ent bands on the Fermi-surface (FS) at sufficiently strong 
magnetic fields. Here we predict a mixing of two or more 
fundamental frequencies in a canonical single-hand Fermi 
liquid with a few extremal FS cross-sections, where the 
MB is non-existent. 

To illustrate the point we consider an anisotropic single 



FIG. 1: Large, 
27rm(^i - 2t)/h 2 . 
Fermi-surface. 



Ai = 2-Km(fi + 2t)/h , and small, A s = 
extremal cross-sections of a layered-metal 



band, with the dispersion, J5(k) , in zero magnetic field, 
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which is a fair approximation for a band in layered metals 
0,0]. Here ku and k± are the in-plane and out-of-plane 
quasi-momenta, t is the inter-plane hopping integral, and 
d is the inter-plane distance. 

When the magnetic field, B, is applied, the spectrum 
Eq.([T]) is quantised as [ijj 

E n (kx) = huj c (n+ 1/2) - 2tcos(kj_d) ±g^ B B/2, (2) 

where lu c — eB cos(0)/m is the cyclotron frequency 
(n = 0,1,2,...), t = t ± J (k F dta,u(Q)) {J {x) is the 
Bessel function), is the angle between the field and 
the normal to the planes, g is the electron g-factor, 
and hb is the Bohr magneton. The s pec trum, Eq.([2]) 
is perfectly 2D at the Yamaji angles [18| found from 
Jo(fc.Fdtan(0)) = 0, where hkp = (2m/x) 1 / 2 is the Fermi 
momentum in pure 2D case, but otherwise there are two 
extremal semiclassical orbits. They give rise to beats 
in dHvA/SdH oscillations with two fundamental FT fre- 
quencies, Fi )S = hAi tS /2ire cos(0) , revealing modula- 
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tions of the cylindrical FS along the perpendicular di- 
rection, Fig.l, as observed e.g. in Sr 2 Ru04 0, [IB|. 

Since there are no different bands one might expect nei- 
ther AB nor MB mixing of the fundamental frequencies, 
Fi and F s in the single-band model, EqJT]) , in contrast 
with canonical multi-band systems 0, Il6j. Actually, as 
we show below, Fi and F s turn out mixed, if n e is con- 
stant, so that a combination frequency F+ = Fi + F s ap- 
pears similar to the AB combination frequency [|[ in two- 
band canonical Fermi-liquids. Using conventional Pois- 
son's summation and integrals [l[ the grand canonical 
potential per unit volume, 
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is its quantum part with the conventional temperature, 
Rt{x) = x/ sinh(x), and Dingle (i.e. collision), < R ^ 
1, damping factors, as derived in Ref. [12]. Differenti- 
ating Q with respect to the magnetic field at constant 
jj,, one obtains the oscillating part of the magnetisation, 
M = -dh/dB, 
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where we neglect small terms of the order of 2t//x <c 
1, and take zero-temperature limit and g = for more 
transparency. 

We are interested in the regime hw c <C 4-7rf, where 
three-dimensional corrections to the spectrum are sig- 
nificant, rather than in the opposite ultra-quantum 
limit [12j, where the quantised spectrum is almost 2D. 
In our intermediate-field regime one can replace the 
Bessel function in Eq. (fl3l) by its asymptotic, Jo (a;) « 
(2/kx) 1 / 2 cos(x — 7r/4) at large x to obtain 
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where B± = Anrnt/ eh cos(0) 3> B. 

Naturally the FT of Eq.© yields two fundamental 
frequencies in the grand-canonical ensemble, where fi is 
fixed, Fig. 2. However, the chemical potential oscillates 
with the magnetic field in the canonical system [l|, ||, 
which affects quantum corrections to magnetisation. Us- 
ing n e = —dQ/dfi, one can find the oscillating compo- 
nent, jl = zlo c /2t: — dfl/d/i, of the chemical potential, 
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FIG. 2: Upper panel: magnetisations as functions of the in- 
verse magnetic field 1/B for grand-canonical (a) and canon- 
ical (b) single-band Fermi-liquids, and their difference (c); 
lower panel: magnetisation FTs for grand-canonical (solid 
line) and canonical (dashed line) single-band Fermi-liquids. 
Here 2nf = w(fi + f a ) = 6500 Tesla, B ± = 300 Tesla, R = 1 
and T — g — 0. 



M = Mo + Aj where = dirh 2 n e /m is its zero-field value 
and 
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Here the "bare" fundamental frequencies, fi tS — m(fj,Q ± 
2t) I eh cos(O) , are now field-independent. Remarkably, 
apart from a normalising factor, the dimcnsionless quan- 
tum correction, z, to the chemical potential, Eq.([7]), turns 
out identical to the magnetisation quantum correction, 
M, Eq.©, which is not the case in a two-band canonical 
Fermi liquid [9(. 
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To get insight regarding the FT of z(B) or M{B), 
Eq.©, we first apply an analytical perturbative ap- 
proach of Refs. [3, LUJ expanding z in powers of R up 
to the second order, z w Z\ + z 2 + z m i X , where 
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yields two first fundamental harmonics with the frequen- 
cies fi and f s identical to those of the grand-canonical 
system, 
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yields two second fundamental harmonics with the fre- 
quencies 2fi and 2f s as in the grand-canonical system, 
and 
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is the mixed harmonic with the frequency F+ = fi + fs , 
which is a specific signature of the canonical ensem- 
ble. Its amplitude is small compared with the first- 
harmonics amplitudes as R(B /2nB±) 1 / 2 in contrast with 
multi-band systems, where the mixed-harmonic ampli- 
tudes have roughly the same order of magnitude as the 
fundamental-harmonic amplitudes (at R = 1) 0, [ill ]. 
Also there is no F_ = fi — f s frequency in the FT spec- 
trum of the single-band canonical system, different from 
the multi-band canonical systems [g, 0] ■ 

To assess an accuracy of the analytical approximation, 
Eg. pop , and some experimental conditions, allowing for 
resolution of the mixed harmonic, we present numerically 
exact magnetisation and their FTs in Fig. 2. Since con- 
vergence of the sum in Eq.0 is poor, one can use its 
integral representation in numerical calculations as 
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where Lt 3/2 (x) = (2/tt 1 / 2 ) dtt 1 / 2 /(e* /x - 1 ) is th e 
polylogarithm. The analytical amplitudes, Eqs. (|8l9|10|) . 
prove to be practically exact with the relative error be- 
low 10 percent at any R, Fig. 3, as the amplitudes of the 
analytical theory of dHvA effect in canonical multi-band 
systems 0, [ll| . Another important feature of the numer- 
ical FT of the solution of Eq. fTTj) is that the resolution 
of the mixed central peak in the middle between two fun- 
damental second harmonics, Fig. 2 (lower panel), essen- 
tially depends on the magnetic-field window used in FT, 
Fig. 4. Since the mixed amplitude is relatively small as 
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FIG. 3: Relative errors of analytical harmonic amplitudes, 
Eqs,(|8l [9l I10[) with respect to numerically exact amplitudes. 
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FIG. 4: Effect of the magnetic field window on the Fourier 
transform of magnetisation. Decreasing the window increases 
the width of FT harmonics obscuring the mixed harmonic. 



(B /2nBj_) 1 ^ 2 *C 1, the window affects its experimental 
resolution. We believe that a relatively small interval of 
the magnetic fields, used in FT, has prevented so far the 
single-band combination frequency to be seen in layered 
metals @, Si, El. Importantly, since the characteris- 
tic field, B± oc Jo(fcFdtan(0))/cos(9), is an oscillating 
function of the tilting angle O, the combination ampli- 
tude also oscillates as the function of the angle, which 
could be instrumental in its experimental identification. 
We notice that the angle dependence of the second funda- 
mental harmonics has been clearly observed in S^RuC^ 

m. 

There is also mixing in the SdH quantum oscillations 
of transverse and longitudinal conductivities. For exam- 



4 



pie the longitudinal conductivity is given by the Kubo 
formula [l9j . 

o = -irhe 2 J dE^p-Tr[S(E - H)v ± S(E - H)v±)], 

(12) 

where v± = 2tdsm(k±d)/histhe longitudinal component 
of the velocity operator, H is the single-particle Hamil- 
tonian including the impurity scattering, and f(E) = 
l/[exp(E — (ij/ksT + 1] is the Fermi-Dirac distribution 
function. Averaging over random impurity distributions 
and approximating the scattering rate by a constant, F, 
one obtains the trace in Eq. (fT2")) as (v±*<s[E — E n (k±) — 
«r] — ) in the ladder approximation. Fhcn applying 
Poisson's summation, one can readily obtain a quantum 
correction, a, to the classical conductivity (for detail see 
Refs. dUMH), which is (at T = 0) 




The asymptotic of the Bessel function, J\{x) « 
(2/itx) 1 / 2 caa{x + tt/4), yields FTs of a very similar to 
those of magnetisation, Fig. 2, with the combination har- 



monic in the canonical system. Generally the scattering 
rate depends on the magnetic field [H| , so that its oscil- 
lations require more thorough analysis of the SdH effect, 
but mixing should be robust. Interestingly, some mix- 
ing of fundamental frequencies may occur even in grand- 
canonical multi- or single-band layered systems, if there 
is an inter-band or inter-extremal cross-section scattering 
by impurities. 

In conclusion, we have found the combination fre- 
quency in the quantum magnetic oscillations of the 
single-band canonical layered Fermi liquid. The differ- 
ence between quantum oscillations of the canonical and 
grand-canonical ensembles, is tiny, Fig. 2, but not ob- 
scured by the MB effect, which is absent in the single- 
band case in contrast with the multi-band systems. We 
have also shown that the analytical (perturbative) FT 
amplitudes are numerically accurate even at zero tem- 
perature and in clean samples (i.e. for R = 1) as they 
are in the multi-band analytical theory [ll[ . A wide 
magnetic-field window is essential for experimental reso- 
lution of the combination dHvA/SdH frequency. 
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